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1 Introduction 

Let {( X ni , Y ni ), 1 < i < n, n > 1} be independent bivariate Gaussian triangular arrays, and let p ni denote 
the correlation coefficient of ( X n i,Y n i ), 1 < i < n. The bivariate maxima M„ are defined componentwise by 



M n — (ALii, AL 12 ) — I max X n j j , max Y n i 


For the asymptotic distribution of M n , Sibuya (1960) showed that M n i and M n 2 are asymptotic independent 
if p n i = p e (—1,1), which coincides with the tail asymptotic independence of Gaussian copula, see Embrechts 
et al. (2002). For the case of p n i = p n , Hiisler and Reiss (1989) derived that 



( 1 . 1 ) 


provided that the following Hrisler-Reiss condition 


lim h 2 n { 1 - p n ) = 2A 2 


( 1 . 2 ) 
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holds with A £ [0,oo] ( the converse assertion is proved by Kabluchko et al. (2009)), where the norming 
constant b n satisfies 

1 -$(&„) = -, (1.3) 

n 

where H\(x,y ), the so-called Hiisler-Reiss max-stable distribution, is given by 

H\{x, y ) = exp ^-4> ^A + e~ v - $ ^A + ^ e~ x ^j 

with $(a;) standing for the standard Gaussian distribution. Obviously, components of M„ are asymptotic 
dependent when A < oo. 

For Hiisler-Reiss model, the asymptotic behavior of the dynamic copula version of normalized M„ has also 
been studied in recent literature. Under the Hiisler-Reiss condition & Frick and Reiss (2013) considered 
the asymptotic behaviors of the distribution of (n(maxi<j< n $(X n j) — 1), n(maxi<j<„ <F(U n i) — 1)). Allowing 
p n i to depend on both i and n, Liao et al. (2014a) extended the result in Frick and Reiss (2013) by assuming 
that 

p ni = 1 - m{i/n)/\ogn (1.4) 

for some positive function m{x). For other work related to Hiisler-Reiss model and its extensions, see, e.g., 
Hashorva (2005, 2006), Hashorva et al. (2012), Hashorva and Weng (2013), Kabluchko (2011), Engclke et 
al. (2014), Das et al. (2014) and reference therein. 

The objective of this paper is to derive the first and the second-order distributional expansions of the 
dynamic Hiisler-Reiss model with p n i given by (11.41) and establish statistical inferences related to the function 
m(x). For the convergence rates and higher-order expansions of univariate extremes, we refer to de Haan 
and Resnick (1996), Nair (1981), Liao et al. (2014b) and reference therein. For the convergence rates of 
bivariate extremes, see de Haan and Peng (1997) for the general case. For the special case of the bivariate 
Hiisler-Reiss model, Hashorva et al. (2014) established the higher-order distributional expansions of M n , and 
Liao and Peng (2014) established the uniform convergence rate of (11.11) . Liao and Peng (2015) also derived 
the second-order expansion of the joint distribution of normalized maximum and minimum. So far, there are 
no studies on the convergence and distributional expansion of M„ under the assumption that (X ni , Y ni )'s 
are not identically distributed. The main goal of this paper is to fill this gap. Borrowing the ideas from Liao 
et al. (2014a), we derive in this paper the limit distribution of the normalized maxima M„ if the function 
m{i/n ) in m satisfies some regular conditions, and establish its second-order distributional expansion 
provided that the convergence rate of maxi< j<„ m(i/n) is given. Furthermore, parametric estimation of 
m{x) is considered through maximum likelihood estimation. The asymptotic properties of the estimators 
can be employed to test the condition proposed by Hiisler and Reiss (1989). 

The rest of this paper is organized as follows. In section [21 we provide the main results and statistical 
procedures. A simulation study and some real data analysis are presented in Section [3] All proofs are given 
in Section [I] 
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2 Methodology 


2.1 Convergence of maxima 


In this section, the limiting distribution and the second-order expansion of distribution of normalized M„ 
are provided with p n i satisfying m- The first result is about the first-order asymptotic which is stated as 
follows. 

Theorem 1. Under the condition Id, 


(i) if maxi<i<„ m(i/n) 0, then for any x,y £ R 

lim P (M n i < b n + x/b n , M n2 < b n + y/b n ) = A(min(x, y)); 

n—>oo 

(ii) if minixi^n m(i/n) —> oo, then for any x,y G R 

lim P (M n 1 < b n + x/bn, Mn 2 < b n + y/bn) = A(x)A(y); 

n—>-oo 

(Hi) if m(s) is a continuous positive function on [0,1], then for any x,y € K 


lim P (M n 1 < bn + x/bn, Mn2 < b n + y/bn) = H(x, y) 

n—> oo 


with 

H(x ’ y)= exp (-«-"/* + iy=|) * -* (vASitj + <*) • 


To establish the second-order distributional expansion of normalized maxima, we consider the follow¬ 
ing three cases in turn: m(t) is monotone and continuous on [0,1]; lim^-nx, maxi<j< n m(i/n) = 0; and 
lim^oo mini <i< n m(i/n) = oo. 


Theorem 2. Under the condition db assume that m(t) is monotone and continuous on [0,1], we have 

lim lQ f n ( P (M n l < bn + x/bn, Mn2 < b n + y/bn) - H(x, y)) 
n ->oo log log n V / 

= \ (l GW)-e <#) v), (2.1) 

where <f(x) is the probability density function of standard Gaussian distribution. 


Theorem 3. Let the norming constant b n be given by m- Assume t/iatlim„_ ) . 00 (logn) 4 maxi< i <„TO(i/n) = 

0, we have 

lim (logn) ( P (M n i < b n + x/b n , M n2 <b n + y/b n ) - A(min(x, y))\ 

= ^ (min(a;, y)) 2 + 2 min(x, y)) e~ min U,y) A(min(:r, y)). (2.2) 

Theorem 4. Let the norming constant b n be given by <d>- Assume that limn^oo (log log n)/( mini<j< n m(i/n)) 
0, we have 


lim (logn) ( P (M n i < b n + x/b n , M n2 < b n + y/b n ) - A(x 

n—>• oo \ 

J2 


2x 


y 2 + 2y 


e v A(x)A(y). 


(2.3) 
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2.2 Parametric inference 


Now we consider statistical inference for fitting a parametric form to the unknown function m(s). Here 
we consider the family m(s) = a + (3s 1 , where a > 0, /3 ^ 0, 7 > 0. Note that when (3 = 0, 7 can not be 
identified, and when 7 = 0, a and (3 cann’t be distinguished, cf. Liao et al. (2014a). 

We use the maximum likelihood estimation (MLE) to get the estimator, which is 


n n I 2 

~ n l0 § 2?r - 2 £ ^ -*&)-£ 2 ( 1 -/) 

2=1 2=1 ' rni/ 


That is. 


(a,P, 7 ) 


is the solution to the following score equations 


In 2 (a, A 7) 
In 3 (<*,/?, 7 ) 


E Tl 

2=1 

E tl 
2=1 



' Pni 

, (1 +pl,)X ni Y ni 

Pndxlt+Y*A 

1 1 

l (logn)(l— 

(log n) (1—p 2 ^) 2 

(logn)(l—P 2 d 2 


( Pni 

, (1 +pl i )X ni Y ni 

P»i( XT+yA) 

1 1 

[ (logn)(l-p* 4 ) 

(log n) (1—p 2 ^) 2 

(logn)(l-p 2 d 2 


[ Pni 

, (l+pli)X ni Y ni 

PniPCb+E, 2 *.) 

1 1 

((l°gri)(l-p 2 i ) 

(log n) (1—P 2 j) 2 

(logra)(l-p 2 i ) 2 


Sr 

i=l 


= 0, 


" Xni * n 


(2.4) 


The following theorem gives the asymptotic normality of the proposed estimator. 

Theorem 5. Assume that m holds with m(s) = a + /3s 7 /or some a > 0, (3 ^ 0, 7 > 0. (ZTien we have 

(2.5) 


A (\/n(d - a), \/n(/3 - /?), Vra(7 - 7)) 4-7V(0,E) 
where the matrices A and E are given by 


(h 


A = 


/( 


0 2(a+/3t 7 ) 2 
1 A 


0 2(a+/3i 7 ) 2 

1 t 7 loff t 


dt 

f 1 t 7 

fl fit 7 log Z 

lO 2(a+/3t 7 ) 2 

JO 2(a+,0Z 7 ) 2 

dt 

r 1 t 27 m 

r 1 fit 2 ^ log £ 

l0 2(a+/3t 7 ) 2 ■ 

lo 2(a+/3tr) 2 

dt 

f 1 i 2 n og i .. 

f 1 /3Z 27 (logZ) 2 

lO 2(a+/3Z 7 ) 2 - 

10 2 (a+j}t*t) 2 


dt 


and 


(Si 


E = 



1 

dt 

r 1 t 7 

dt 

r 1 £ 7 log £ 

rl+\ 

2( 

Q! + ,flZ 7 ) 2 

JO 2(c*+/3Z 7 ) 2 

JO 2(a+fit-r)' 

jLLi | 


t 7 

dt 

r 1 t 27 

dt 

r 1 £ 27 log £ 

-ri+ 

2( 

a +/3Z 7 ) 2 

lo 2(a+,8Z 7 ) 2 

0 2(a+/3Z 7 V 

jtLL 


£ 7 log £ 

dt 

r 1 Z 27 logZ 

dt 

fl C 7 (logt)^ 


2( 

ct+,i3z 7 ) 2 

JO 2(a+fit 7 ) 2 

JO 2(a+Pt 7 )' 



( 2 . 6 ) 


Another interesting parametric form is m(s) = a + /3s for some a > 0, (3 £ R. In this case, when 
/3 = 0, m(s) becomes constant, which means that the observations - ■ ■ ,(X n ,Y n ) are independent 

and identically distributed random vectors. 


Theorem 6 . Suppose (11.41) holds with m(s) = a + /3s for some a > 0, (3 ^ 0. Then we have 

( 0 , e) , 


^ (A^ (‘ + i) - (ifSSs + w >'08 (! +1))) 


JfiiX _ 4- log (l + 4") - (-4- + -ASA: 

A \2/3 2/3 2 S V “/ \2/3 2 ^ 2/3 2 (a+ 

where E zs given by 


— pot 

+P) 


2 &§— 

2/9 : 


^i<.g(i+4))), 


A 


(2.7) 


E = 


2o(q+/3) 

~ 2/3(cl+/3) + W 7 l0g ( 1 + a) gk( 1 + 


1 + # 


2/3(i+/3) + W 1 1°S 

log (l + |)) 
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3 Simulation and data analysis 


In this section we examine the finite sample performance of the proposed estimators by drawing indepen¬ 
dent (X n i, Y n i), • • • , ( X nn , Y nn ) with (X n i,Y n i) following the bivariate Gaussian distribution with coefficient 
p ni = 1 — m(i/n)/\ogn. We consider n = 1000,3000 or 10000, and repeat 1000 times. 

First we consider m(s) = a with a = 1 or 10, and calculate the average and mean squared error for a. 
We can observe from Table [T] that i) the averages of a is near by the true value a; ii) small mean squared 
errors show the robustness of a. Next the case of m(s) = a + /3 s is considered. Table [5] reports the averages 
and mean squared errors for estimator (&,$)■ As n becomes large, the accuracy of all estimators improve. 
Finally, we consider the case of m(s) = a + /3s 7 with sample size n = 10000. The simulation shows that all 
estimators are closer to their true values with small mean squared errors, cf. Table [3] for details. 


Table 1: Estimators for the case of ro(s) = a 



a = 1 

n = 1000 

a = 10 

n = 1000 

a = 1 

n = 3000 

a = 10 

n = 3000 

a = 1 

n = 10000 

a = 10 

n = 10000 

E(a) 

MSE(a) 

0.9980043 

0.002109944 

9.992257 

0.02555685 

1.001875 

0.0006846334 

9.997267 

0.01812748 

1.00005 

0.0002033414 

9.999436 

0.007965965 



Table 2: 

Estimators for the case of 

m(s) = a + /3s 

with a = 1 . 



P = l 

o 

II 

P=1 

o 

II 

(3 = 1 

o 

II 

CCL 


n = 1000 

n = 1000 

n = 3000 

n = 3000 

n = 10000 

n = 10000 

E(d) 

1.002458 

1.002368 

1.001585 

0.9973226 

0.9996816 

1.000073 

MSE(d) 

0.01125451 

0.00864667 

0.00385037 

0.002605236 

0.00124363 

0.000797927 

E(/3) 

0.9978728 

-0.000328 

0.9967 

0.004754433 

1.001072 

0.000712807 

MSE(/3) 

0.04707124 

0.02588997 

0.01669121 

0.007914633 

0.00509091 

0.002447169 


Table 3: Estimators for the case of m(s) = a + /?s 7 with a = (3 = 1. 



E(d) 

MSE(d) 

E(/3) 

MSE(/3) 

E(7) 

MSE(7) 

7 = 0.5 

0.9437925 

0.05247714 

1.058977 

0.04504459 

0.5092629 

0.03055645 

7 = 1 

0.994503 

0.004061697 

1.009859 

0.005172654 

1.01609 

0.03928909 

7 = 1.5 

0.9955838 

0.001943133 

1.003917 

0.005124053 

1.504647 

0.06171575 


For the applications, we consider four couples of real data sets: The first is the log-returns of the exchange 
rates between US dollar and British pound and those between Canadian dollar and British pound from April 
3, 2000 to November 11, 2014. The second is the log-returns of the Shanghai Stock Exchange composite 
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index (SSE Composite) and ShenZhen Stock Exchange Composite index (SZSE Composite) from March 4, 
1996 to November 12, 2014. The third is the log-returns of the CSI 300 index and CSI 300 index futures 
from April 16, 2010 to November 13, 2014. The forth is the wave and surge heights in southwest England 
which comprise 2894 wave heights and 2894 surge heights. All time series are plotted in Figure Q] 

First, we calculate the *th sample correlation for each couple of the mentioned data sets by using 
{(X n i, Y n i), ■ ■ ■ , (X n i,Y n i)}. Figures EE] show respectively that each tends to constant ultimately. Now 
we estimate the correlation p = 1 — m(i/n)/ log(n) by assuming that m(s) is a constant, which also are 
illustrated by Figures EMI respectively. The constancy of m(s) shows that observations are identically 
distributed. 


Log returns 


Log returns 



Log returns 



Wave and surge 



CSI300 index 


Wave 


Figure 1: log-returns of the exchange rates between US dollar and British pound and those between Canadian 
dollar and British pound (top left); the log-returns of SSE Composite and SZSE Composite from (top right); 
the log-returns of the CSI 300 index and CSI 300 index futures (bottom left); the wave and surge heights in 
southwest England (bottom right). 


4 Proofs 

The aim of this section is to prove our main results. In the sequel, let Fi(x,y) denote the distribution 
function of ( X n i,Y n i ), 1 < i < n; and let u n (x) = b n +x/b n for notational simplicity. 

Proof of Theorem [l] We only consider the case (iii) here, since the other two cases can be derived by 
Slepian’s Lemma and the result of case (iii). 
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Figure 2: Exchange rates. Dotted line represents the sample correlations, and solid line represents the 
correlation estimate p = 0.4738478 with m = 4.338189. 



Figure 3: SSE Composite and SZSE Composite. Dotted line represents the sample correlations, and solid 
line represents the correlation estimate p = 0.9048648 with rh = 0.8009351. 
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Figure 4: CSI 300 index and CSI 300 index futures. Dotted line represents the sample correlations, and solid 
line represents the correlation estimate p = 0.9455578 with m = 0.3817058. 



Figure 5: Wave and surge heights. Dotted line represents the sample correlations, and solid line represents 
the correlation estimate p = 0.2955482 with to = 5.614759. 








It follows from (11.31) that 


, /01 a log log n + log 4n 

b n = (2 log n) =-:-—-h o 


2(2 log n) i V(logn)^ 

which implies that 6 2 ~ 2 logn as n —> oo. Combining with (11.41) . we have 
U n (x) - p n iU n (z) 


Pni 


1 Pni 

1 + Pni ' b n ^J\ - f? m V 1 + Pnib: 


X — Z 


1 Pni Z 


m(-) 

v n / 


\ ( lo g™) ( 2 - Sir) 


1 log log n + log 47 r 

(2 log n ) 2 1 - - 2 —-—+o 

' 4 log n 


logn 


(4.1) 


( 2 logn )2 ^1 — 


log log n+log 47r 
4 log n 


( 1 )) /^c 

\ log n ) j W (log n) y log n J 


m(-) 

' n. / 


4 - _hl:___ 

\ (log n) ( 2 - ^p) (2 logn)^ (l 


log log n+log 47T 
4 log n 


V log n J J 





x — z log log n f i \ log log n x — z 


2 ^/m (-) 41og n\] n \n) + 81ogn ^/ m (i.) + ° ^logn J + U ^logn 


+ 0 


( Z 


(4.2) 


for large n. 

By using the inequality |<t>(:r) — c E > (j/)| < \x — y\ for any x, y € M, we have 

/ 


$ 


U n {x) - pniU n {z ) 


V 1 ~ P 


- $ 


\\ 


m ( - ) + 

n 


2 \/ m (r) 


<(i + N)o 


log log n 
log n 


+ (i + N)o 


logn 


for large n and any z £ R, which implies that 


exp —z - 



2 bl 


dz 


2 \J m (n) , 


e Z dz{ 1 + o(l)) 




x — z 


e Z dzdt 


2 s/ rn (t) / 

= -e~ x + e-y I $( v / ^)+ X ZdL ) dt + e~ x [ f v'mjtj + V ~ ' r 


as n —> 00 . 

Note that by Castro (1987) and (11.31) . 


2 y/m(t) 


tp(b n ) 


2y/m(t) 


dt (4.3) 


n-i = l-$(6„) = ^(l-6- 2 + 0(6- 4 )) 

On. 


for large n, and Nair (1981) showed that 


lim bl (—n( 1 - $(u„(x))) + e x ) = X S~ e 1 

n—>00 Z 


(4.4) 


(4.5) 
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Combining with (14.31) . we have 

n 

-^2 (! - Fi(u n (x),u n (y ))) 

2=1 

n poo / 

■ (1 - $(m„(x))) - n~ l (1 + b~ 2 + 0{b ~ 4 )) J2 $ 

2=1 J y \ 


= —n (1 — 


U n (x) - PniU n {s)\ ( S* . 

exp —s — I ds 




Pli 


262 




S dsdt 


2 v / m(t ) } 

= -e-® / $ ( + df-e“ x /% + y_:r 


J o \ 2yjm(t) j 

as n —> oo, which implies the desired result. 
The proof is complete. 

Proof of Theorem [2j By (14.21) we can get 


2 \Jm{t) 


dt 


(4.6) 


□ 



e 


/ 1 


log log n 


\logn 


41ogn 
log log n _ x f 1 


L l K' A ^ + ^fe)( v ^' 


x — z 


2 yjm{t) j 
Vm(t)ip ( yjm(t ) + y 


2y/m(t) 


e Z dzdt 


dt 


2 log n J o y 2yjm(t) / 

as n —> oo. 

By Taylor expansion with Lagrange reminder term, we have 
M n (x) - p ni u n (z) \ 


(4.7) 


$ 


V 1 - p 






x — z 


2 \l m ( j k). 


where 


U n (x) - p ni U n (z) I ( i\ x-z \ I U n (x) - PniUn(z) / / * \ X~Z 

mm | - \ m — M--== < 6 < max —-— —>/m - + 


V 1 ^ /4 

Combining with (14.71) we have 


2 V m (£). 


a / 1 “ P 


2 V m (£>, 


lv^ f°° I ~ ( U n (x) - PniU n (z ) \ „ 1 \ X - z 




V 1 - pi 



— $ | \ / m ( — ) + _ 

nJ 2 7^y, 


e 2 dz 
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Jy 


1 n r°° 

~7tYJ M <0 


2 V m (-), 

u n (x) - p ni u n (z ) 


U n (x) - PnjU n (z) 

V 1 ~Pli 



— \ m 


x — z 




Z dz 




log log n ^ 


°ni 


2 log n 

as n —> oo since 


-1 rL pOO 


V 1 - p 

Vm(t)ip ( \Jm(t) + V 



x — z 




e Z dz 


2 y/m(t) 


dt 


(4.8) 


om 


U n {x) - p n jU n (z) 

V 1 ~Pli 



— \lm 


x — z 




e~ z dz = O 


( / loglog i 

l V i°s n 


' »=i J y 


Now, we first assert that 

-4 e~ z dz - 

2 V TO (n)/ 

holds for any x, y £ R. Combining with (14.81) . we can get 

1 V- / 



<J> yjm(t ) + 


0 Jy 


2yjm(t) 


e Z dzdt = O ( — ) (4.9) 


— r °° ( U„(x) - pniUn(z) 


hJy 


3 li 


V V 1 - p 


e Z dz + 


$ y/m(t ) + 


x — z 


y-x 


2 y/m{t) 


0 Jy 


dt 


2y/m(t) 


e Z dzdt 


(4.10) 


as n — > oo. 

From (14.41) . (14.51) and (14.101) . it follows that 


n pi poo / 

-V](f - Fi(u n (x),u n (y))) + e~ x + / 4> yj m(t) + 

i=i •'o •'a V 


Z dzdt 


n p OO 

= —n(l - $(«„(*))) + e -x - n" 1 (l + 6“ 2 + O (6" 4 )) £ / 4> 

i=i J v 


2\/m(t)) 

U n {x) - PniUn(z) 


V 1 ~ P 


e z exp l — 


z 

26S 


dz 




x — z 


e Z dzdt 


2 yjm{t) 

-j n p OO / 

= -n(l - $(u B (x))) + e _x -4> 

n i =l ^ V 


Mn(^) ~ PniUnjz) 

V 1 - pL 


* dz + Jo J™ + 


2 \W) 


nh l ht-ly 

log logn^_ x 


$ 


- PniUnjz) 
K \A-Pni 


1 — ) dz + O (6„ 4 ) 


21 ogn 

as n —» oo, which implies that 


\Jm[t)ip yj m[t) + 


y - x 


2 yjm(t) 


dt 


P(M n 1 < u n (x),M n 2 < u„(y)) - H(x,y) 


Z dzdt 
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= H{x,y) ^exp ^2\ogFi(u n (x),u n (y)) +e x +J J $ (^/m{t) + 

( n pi pOO / 

~E (1 - Fi(u n (a;), u„(y))) + e~ x + J J $ ( yjm{t) + 


x — z 


2 \A n{t) 


e Z dzdt 


x — z 


2y/m(t) 


- 9 (! - Fi{u n (x), u n {y))f (1 + o(l)) J (1 + o(l)) 


2=1 


Ti Sr(i 


y-x 


2 V m (t) 


dt ] e x H(x, y) 


as n —> oo. 


The remainder is to show that holds for any fixed x,y G R. Without loss of generality, we 
that m(t) is increasing. 

If x < y, note that ^ \/m(t) + ^ e~ z dz is increasing about t, so we have 


and 


1 n pOO ( 

=£/.*( 

tjtj 

J —— J V 


x — z 


2 \) m ( n ) , 


e Z dz 


$ I \lm ( - 

n 


x — z 




Z dzdt 


<*— pOO 


> EL / $ + 


2=1 


2 y / m(t) 


e Z dzdt 


fjy(vw) + 


2 


Z dzdt 


E/ * 


n 


2=1 17 V 

i+1 


x — z 




e 2 dz 


= E 


i=i J ; 
r-l + i /-oo 


$ I \lm [ - 

n 


x — z 


2 \J m (n) , 


e Z dzdt 


< 


j '[% «+ 


'0 •/!/ 
f*l /*oo 


2 \/m(t) 


e Z dzdt 


o 


$ ( \Jm{t) H- ’ ] e 2 dzdi + O (— 


which implies that (Id.91) holds for a; < y. 

To verify (|4.9|) holding for a; > y, we just need to prove that 

r 1 rv 


1 


E * 


I m 


2=1 


— J j e 2 dz — / / $ f i/m(t) H '* I e Z dzdt = O 

2y^iy/ Jo jX V 2 ^/m(t) / \n 



assume 


, (4.11) 
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which will be proved in return by the following three cases: (i) y < x — 2 m(l); (ii) x — 2 m(l) < y < x — 2 m( 0 ), 
and (iii) x — 2m(0)<y < x. In fact, the arguments of (i) and (iii) are similar. The rest is to focus on (i) and 
(ii). 

For case (i), i.e. y < x — 2m(l), it is known that y < x — 2m(t) < x for any t £ [0,1]. Hence, 


and 


1 " r x ~ 2m («) 

71 < ^ 


i =l J v 


$ I ( /m| - 
n 




e Z dz 


- e/: 

i— 1 J ~n~ 

n „i_ 

< zjL 

_ n _ z*^: /* 

< E/„/ 




$ Wm 


x — z 


2 \/ m (i), 


Z dzdt 




$ y/rn{t) + 


2 \/m(t) 


Z dzdt 


x—2 m{t) 


2=1 

1 fX—2 m{t) 


•h ( \/m (f) H- X ) e Z dzdt 

\ 2^w) 


0 Jy 


$ ( \Jm (t) H- 2 ) e Z dzdt 


71 ‘ 


*=1 


$ | (/m|- 
n 


x — z 


2 vMf), 


> 


> 


zjrf 

i=i J n J v 

n r ^ r 

E/ 


i±i ra-2m(i) 


$ 


X — Z 


2 ]/ m (i), 


iti i-x-2m(t) 


$ \ m 


x — 0 


2 \Mf), 


e Z dzdt 


e Z dzdt 


i±i r x-2m(t) 


E/ 7 / 

. — 1 */ — -'l/ 


1 nx—2 m(t) 


$ ( \/ m (i) H- ' r . ] e 

V 2^/rnJt) J 


0 Jy 


$ ( \/m (f) H- r ] e Z dzdt + O 

y 2v^rm/ 


(4.12) 


(4.13) 


Similarly, 


71 < ^ 


_j ./ x — 2m(^ ) 


and 


$ \ m - 


2 V to (f), 


2 dz < / / $ ( v /m (t) H- * | e z dzdt+0 ( — 

Jo Jx— 2 m(t) \ 2y/m (t) J \n 

(4.14) 


i E 

71 ' 


$ ( m - 


*=i ^ x - 2 m (4:) 


1 pX 


2 \/ TO (n), 


e 2 dz > 


0 J x—2m(t) 


$ \/m (f) + 


x — z 


2 \Jm {t) 


e Z dzdt.( 4.15) 


Combining with (14.1211 - (14.151) . it shows that (14.111) holds for case (i). 


13 


















































Next we consider case (ii), i.e. x — 2m(l) < y < x — 2m(0). Note that there exists x* £ (0,1) such that 
y = x — 2 m(x*) since m{t) is increasing and continuous. Split the following integral into two parts: 


pi px—2 m(t) / 

Jo J i $ ^ + 


x — z 


2 \Jm (t) 


e Z dzdt 


x* i'X—2m{t) 


$ \Jm (t) + 


x — z 


/o jy \ 2sjm (t) J 

By arguments similar with (14.1211 - (Id. 151) . we can get 


dzdt - 


pi px—2m(t) / 

/ / $ i/m(t) + 

-'z* Jy \ 


X — Z 


2 \/m (t) 


Z dzdt. 


IV l./m" 

71 . 


1=1 J y 

px* px— 2 m{t) 

o *4 


x — z 




e Z dz 


and 


1 E 

n < ^ 


<f> (f) + 

- 2 m(^ ) 


- = ) + o(- 

2y/m (t) J \n 


i=[nx*]-\- 1 1 


$ I Wm| - 

n 


x — z 




Z dz 


pi px—2m{t) / 

J J $ Em (t) + 


' = ) e" 2 cfadt + O (- 

2yjm (f) J \n 

Combining above with (14.141) . (14.151) . we show that (14.111) holds for case (ii). 

Now, (14.111) is derived for any fixed which complete the proof. 

Proof of Theorem [3j For fixed x, y £ R, if max(j:, y) < z < 4 log b n we have 


□ 


^ / Un(min(s,;/)) - Pni'U'nfe') | ^ 


exp - 


(_&i 


— pni) _ min(x,^) p ni z 


1+Pni 1+P: 


V 1 ~ P 




for large n by using Mills’ inequality. Combining with ED , ED and Hindoo (log n ) 4 maxi<i<„ m{i/n) = 0 , 
we have 

f 4 '°“’' 4 ( e -exp (-|1 | dz 


K(x,y) 


< 


< 


(1 +Pni) exp 


—pni) min(x, 2 /)+max(ai 

T+JTn 


,y) j 


2\/2m (i) (l - loglos "+ log47r 


fK~ ( max(q,,?/) —min(x,p) _ / l-p n i _ l l-pni 4 log j 

V l V 1 ^' " V^i ^ 

(logn)) eX P (E TO (£> ( 


4 log n 


log log n+log 4tv 
2 log n 


.p_Y) 

V log ™ / / 


+ 


\x+y\ 
o Mi/n) 

Z ~ logn 


max(x, y ) - min(x, y) - 2 m(i) (l 


log log n+log 47 t 


2 log n 


+ O 


(^J\) 

V log n J J 




i\ ( log log n + log 47 t / 1 

< 2\I2 max m — 1 -—- 1 - o - - 

1 <i<n \nJ \ 4log n \log?z 


exp - 5 mmi<i<„ m 


(i)( 


i _ log log n+log 4 tt . 

1 2 log n ' U 


(^J\) 

V log n J J 


+ 


\x+y | 


2 -- 




log r 


(max(i, y) - min(x, y) - 2 m (+ ( 


i _ log log n+log 4 tt i 
' 2 log n ^ 


V log n J J 


_ 4 maxK^ m(,yn) bg ^ 
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= 0 ( 6 - 4 ) 

for any 1 < i < n. 
Noting that 


/; .( 

J 4 log b n y 


Mn(min(x,y)) - p n jU n {z) 

\A - Pni 


e 2 ex P ( -^2 ) dz = °(0 


for 1 < i < n, we have 

n 


$ 


2=1 


jc(x,y) 


u n (min(a;,y)) - Pni'U'n (z) 

V' 1 ~Pni 


e exp — 


2 b 2 


dz = 0{b n ) 


for large n. Hence combining above with (14.51) . we can get 
b n (- £ (1 - Fi(u n (x),u n (y))) +e-““(x,v)^ 

= b 2 n (-n(l - $(min(x, y))) + e~ min ( x ^ 


n 

- bln - 1 (1 - 6" 2 + 0(6- 4 )) -1 / 

i= l 

(min(x, y)) 2 +2 min(:r, y) _ minfx ,^ 
2 




m(x,v) 


as n —> oo, which implies (12.21) . The proof is complete. 
Proof of Theorem [4J By Mills’ inequality we have 


u n (mm(x,y)) - p n jU n (z ) 

V V 1 - P 2 ni 


e z exp i — 


z 

2ft 


1 - $ 


U n (y) - PniU n (z) 


V 1 ~ P‘ 

for large n , which implies that 


< 


eX P (- Vl+ftTo “ J l+iSf “ 5 - Pni)) 


^( 1 + ^ ( i4„o + <t) 


r41 ° Sbn ( ! _ $ [“»(*0 - 


's/I - P 


e z exp i — 


26S 


dz = 0 ( 6 " 4 ) 


for any 1 < i < n due to (fT~H) . (IXT1) . lim^oo mini<i< n to (i/n) = oo and lim^oo = C 


Combining with 


1 - $ 


/ 4 log b n 


Un(y) - PniUnjz ) 


e exp 


26S 


dz = 0(b n ) 


we have 


= n 


a 

^ ^ 'll*nipt)^ni ^ ^n(2/)) 

n p oo / 

(l- 6,- 2 + 0 ( 6 - 4 ))- 1 ^ / 

*=i V 


i=l 

-1 


1 - $ 


Mn(i/) - PniU n {z) 

, V 1 - Pni , 


e exp 


z 

2 K 


dz 


= 0 ( 0 - 


dz 


□ 


(4.16) 
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It follows from (14.51) and (14.161) that 


M F i{u n {x),u n {y)))+e x + e y 


i= 1 


= b 2 n ( -n(l - $(u„(x))) + e x - n(l - $(u„(y))) + e v + ^P(X ni > u„(x), > u„(y)) 


i=l 


x 2 + 2 x 2/ 2 + 2 y 

2 2 

as n —► oo. Hence (12.31) can be derived, which complete the proof. 

Proof of Theorem [5j Define 

1 + Pni 


Zi = — 


(log n) ( 1 -/ 4 )' 
“ Zi t i + Zj,2 + ■Z’i,3) 


(-^"ni + Yni) + 


(logn)(l — p„j) 


2 \2 XniYni 


(log n)(l - p 


one can check that 


E^i 

E^ ? 3 

e z f _ i Z: L 3 


_ 4(p^+2p 2 ,) 

(logn) 2 (l-Pni) 4 ’ 

_ _ Pm 

(log n) 2 (1—.) 2 ’ 

2 Pn 


EZ 2 


— _ ( 1 +pL) 2 _ /I _i_ 0 0 2 \ 

floe - n'l 2 T1 — o 2 .'I 4 ' rni) 


E ZipZip = — 
E /(, 2 a = 


iill 

(log n) 2 (l-p 2 
6pL.(l+P 2 ,) 


(logn) 2 (l-p 2 i ) 4 

... _ pL.(i+pD 

(logn) 2 (l—p 2 ^) 3 ’ i,3 (log n) 2 (l — p 2 ^) 3 ’ 


which implies that 


-j n i n 

-Ve^ = iv 

n ^ r? ^ 


i=l 


i=l 


/ i °+K*r i (°+*(*r ) 3 ^ 

log n ' 2(log n) 2 

( 2 (“ + ^(h 1 ) ! (i- f ^ : )') 


/o 2 (a + /if 7 ) 2 




as n —> oo. It is easy to check that 

E (y E ( z - - E z ?)) = ^ X> ( z . 2 - E Z ) 2 = ° Q). 

which combining with (14.171) implies that 

Obviously, we have 


max 


1 


\/n 


l<i<n 

Hence combining with (14.18[) . (|4.19l) we can get 

1 


0 and E max — Z~ = o( 1). 

1 l<i<n Tl 7 


2(a + W 


1 


dt 


□ 


(4.17) 


(4.18) 


(4.19) 


(4.20) 


as n —> oo. 
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Define 


2 ? = — 


{^Pr* 


(log n) (1 — PniY 


+Yni) + 


i^y^+pii) x y («) 7pni 

-A-n.i 1 77/ 


(logn)(l — p: 


,2 ^2 


:= z; 


2,1 


J i,2 


■z: 


7,3’ 


and 


zr = - 


(^) 7 ( 1o S^)p™ ^2. + y2J + (^) 7 ( lQ g^) ( 1 + pli) 


(log n) (1 Pni) 

_ ry*- t= . . ry** 

1 "i" ^7,2 “i" ^7,3* 

Similar to the proofs of (|4.20|) , we can show that 


(logn)(l-p2.)2 


(logn)(l-p2.) 


(^) 7 ( lQ g^)Pm 

(logn)(l -p2.) 


lim„^ TO l n2 {a, (3, 7 ) — fo, 

linin-z-oo y^}n 3 ( 01 , J3,j) = N ^0, f 0 2 ( a +pry)* 


By arguments similar to (14.171) . we have 


1 n 1 n 

lim -VEZiZ;= lim -V (-) E Zf = 

n —^00 Tl z ' n—too 77. ^' \ 77 / 


7=1 


7=1 






-76 — 76 

lim -VEZ*Z** = lim -V 

^^ n—>oo 77 , ^' 


n—>oo fi 


7=1 


7=1 


and 


-V EZiZ?* = lim - V 

n—>■ oo Tl ^^ n—>■ oo Tl z ' 

7=1 7=1 


lim 


2 7 


log- 


0 2(a + /3t 7 ) 2 

rl t 27 log t 


EZf = 


2(a + /3f 7 ) 


rdt 


log; 


Hence, by Cramer device, we can derive that 

1 


E Z~ = 


lim ~^= (l n i(a, /3, 7 ), Z n 2 (a, 13, 7 ), Z„ 3 (a, A7)) T = AT (0, E). 

i->oo ^/ri 


i 7 log t 
2(a + /3t 7 ) 


rdt. 


where E is given by EH). 

It is straight forward to check that 


lim — 

n — ^00 Tl 


dlnl(<*,Pn) din l(Q,ft,7) 


a 1 

£ 


da 

9/3 


d-y 


C^n2(c»=,/3,7) 

dln 2 (aj 3 , 7 ) 

dl 

r.2(a,/3,7) 


da 

9/3 


d-y 


din 3 (a,/3,7) 

a/„ 3 (a,/ 3 , 7 ) 

dl 

n3(“./3,7) 


\ da 

a/3 


07 / 


1 rtf 

r 1 * 7 

dt 

pi jSi'Mogi 

2(a+/3t7)2 UI J 

0 2(ct+/3tr)2 

Jo 2 (ct+atr ) 2 

* 7 rjf 

r 1 i 27 

dt 

pl /3lfc 27 log£ 

2(a+^t7)2 al J 

0 2(a+/3f>)2 

Jo 2(a+at7 ) 2 

* 7 log * ,7+ 

rl t 27 lo S Z 

dt 

fl /3t 27 (logt ) 2 

2(a+/3t7)2 u '' J 

0 2(ct+/3fT ) 2 

Jo 2(a+at7)2 


:= A. 


(4.21) 


(4.22) 


(4.23) 


Hence, the desired result is derived by (14.221) . (14.231) and Taylor expansion. The proof is complete. □ 

Proof of Theorem [6j It follows from the proof of Theorem [5] with known 7 = 1 . □ 
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